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Abstract 



We investigate the disappearance of collective flow in the reaction 
plane in heavy-ion collisions within a microscopic model (QMD). 
A systematic study of the impact parameter dependence is performed 
for the system Ca+Ca. The balance energy strongly increases with im- 
pact parameter. Momentum dependent interactions reduce the balance 
energies for intermediate impact parameters b ~ 4.5 fm. Dynamical 
negative flow is not visible in the laboratory frame but does exist in the 
contact frame for the heavy system Au+Au. For semi-peripheral colli- 
sions of Ca+Ca with b ~ 6.5 fm a new two-component flow is discussed. 
Azimuthal distributions exhibit strong collectiv flow signals, even at the 
balance energy. 
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I. INTRODUCTION 



The prediction of collective flow in heavy-ion collisions by the hydrodynamical 
model has yielded a powerful tool for the investigation of excited nuclear matter. 
Main goals are to determine the equation of state (eos) and the in- medium nucleon- 
nucleon cross section. One possible approach is the measurement and calculation of 
the transverse flow in the reaction plane. At beam energies above 100 — 200 AMeV 
two-body collisions rule the dynamics yielding the typical bounce-off behaviour 0-0, 
which is the deflection of cold spectator matter from hot compressed participant mat- 
ter. The attractive part of the mean field becomes more and more important with a 
decrease in energy. As a consequence even negative scattering angles are possible || 
which can be imagined as partial orbiting of the two nuclei 0]. At a certain incident 
energy, called the balance energy E ba i, the attractive and repulsive forces which are 
responsible for the transverse flow in the reaction plane cancel each other, causing 
the disappearance of this particular flow characteristic. 

The notation "energy of vanishing flow", as the balance energy is often called, can 
lead to misunderstandings: In particular, we will demonstrate by inspecting az- 
imuthal distributions that strong flow still exists at the balance energy. Whereas it 
was shown for small impact parameters that the balance energy depends only weakly 
on the stiffness of the equation of state [§|§, a large sensitivity to the nucleon- 
nucleon in-medium cross section was recognized H|J|7[]. The functional dependence 

of the balance energy on the system size can be approximately described by a power 

_ i 

law: Eb a i ~ A t J |T0|j7||. Systematic studies of the mass dependence of the disap- 
pearance of flow proposed a reduction of the in-medium cross section of about 20% 
with respect to the free NN-cross section at normal nuclear density [7J] by comparing 
the measured data |Tl]-|T6|,|7| with BUU calculations. However, all investigations ne- 
glected to study the impact parameter dependence of the disappearance of flow. 
In this contribution we show that a variation of the impact parameter changes deci- 
sively the balance energy E ba i(b) and as a consequence the mass dependence analysis 
receives an important new variable. 

The system Au+Au exhibits no negative flow in the laboratory frame. However, if 
the initial pre-contact rotation of the system due to Rutherford-trajectories is sub- 
tracted, large negative flow appears. 
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A new two-component flow appears in collisions with large impact parameters. 
Azimuthal asymmetries persist at the balance energy. 

The balance energy E^i is nearly independent of particle type @], although it is well 
known that the strength of the flow depends on it. Therefore we will mostly regard 
all nucleons and check the effect of taking clustering into account. 



II. THE MODEL 

The Quantum Molecular Dynamics model (QMD) [pi[l7H2^| is employed here. 
In the QMD model the nucleons are represented by Gaussian shaped density dis- 
tributions. They are initialized in a sphere of a radius R = 1.12A 1 / 3 fm, according 
to the liquid drop model. Each nucleon is supposed to occupy a volume of h 3 , so 
that the phase space is uniformly filled. The initial momenta are randomly chosen 
between zero and the local Thomas- Fermi- momentum. The Ap and At nucleons 
interact via two- and three-body Skyrme forces, a Yukawa potential, momentum 
dependent interactions, a symmetry potential (to achieve a correct distribution of 
protons and neutrons in the nucleus) and explicit Coulomb forces between the Zp 
and Zt protons. Using this ansatz we have chosen a hard equation of state with a 



compressibility of k = 380 MeV p3| , |2^| . For the momentum dependent interaction 
we use a phenomenological ansatz which fits experimental measurements 

28|,|27j of the real part of the nucleon optical potential. The nucleons are propagated 



according to Hamiltons equations of motion. A clear distinction is made between 
protons and neutrons with Coulomb forces acting only on the protons and an asym- 
metry potential containing the asymmetry term from the Bethe-Weizsacker formula 
acting between protons and neutrons. Furthermore parameterized energy dependent 
free pn and pp cross sections are used instead of an averaged nucleon-nucleon cross 
section. They differ by 50% at 150 MeV. It was shown that their energy dependence 
cannot be neglected |[29|| . Hard N-N-collisions are included by employing the collision 
term of the well known VUU/BUU equation P, p4| , p0| -|33|] . The collisions are done 



stochastically, in a similar way as in the CASCADE models |34] , |3"5|j . In addition, the 
Pauli blocking (for the final state) is taken into account by regarding the phase space 
densities in the final states of a two body collision. 
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III. RESULTS AND DISCUSSION 



For the investigation of transverse flow in the reaction-plane the in-plane trans- 
verse momentum p x is usually plotted versus the normalized rapidity y/y p . Fig.l 
shows the p x {y) distribution at two different energies for the system Ca+Ca and 
b = 0.5b max « 4 fm. At 80 AMeV a negative slope (corresponding to negative scat- 
tering angles) is visible whereas for 130 AMeV the opposite sign (positive scattering 
angles) is found. The first corresponds to negative scattering angles of the majority 
of the protons, the latter illustrates the deflection of nucleons caused by nucleon- 
nucleon collisions. 

In order to determine the balance energy, the energy is varied between these two 
values and a linear fit is applied to the slopes of the p x (y) distributions. These 
slopes, which are called reduced flow, have negative values for energies smaller than 
Eb a i and positive values for energies higher than E^ai- The balance energy E^ai 
is obtained again by a linear fit to the energy dependence of the reduced flow at 
the point where the reduced flow passes through zero (fig. 2). Onethousand events 
of Ca+Ca are performed for a hard equation of state without momentum depen- 
dent interactions. Different symbols correspond to the different impact parameters 
0.25b max , 0Ab max , 0.5b max , 0.6b max . The balance energies differ completely for the dif- 
ferent impact parameters. This is in contrast to claims in p6| . The errors of the 
balance energies are approximately ±5 AMeV . 

Fig. 3 depicts the impact parameter dependence of the balance energy for the system 
Ca+Ca. An approximate linear increase of the balance energy with impact parameter 
is visible. At larger impact parameters fewer nucleon-nucleon collisions yield reduced 
repulsive forces, therefore the attractive mean field dominates. For larger impact 
parameters the balance energy is smaller if momentum dependent interactions (mdi) 
are included, due to their repulsive effects. The balance energy is insensitive to the 
inclusion of mdi for small impact parameters b < 0.25b max . The balance energy for 
Ca+Ca varies from 65 to 150 AMeV without mdi and from 75 to 115 AMeV with 
mdi, depending on impact parameter. Experiments show the balance energy for 
Ar+Sc , i.e. A = 85, to be 87 ± 12 AMeV, the impact parameter was estimated 
to be approximately 0Ab max ~ 3fm. This value is compatible with ours. Even for 
rather central collisions with a maximum impact parameter of 0Ab max the balance 
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energies for Ca+Ca reach values from 65 AMeV up to 95 AMeV depending on impact 
parameter. This is a significant variation contrary to the claims in |]3l)fl . A precise 
knowledge of the impact parameter is of utmost importance before any conclusions 
about the balance energy concern the equation of state or the in-medium nucleon- 
nucleon cross section. 

Let us now turn to a different question: Two flow-components appear in one event 
showing both positive and negative flow if semi-peripheral collisions of Ca+Ca at 
b = 0.85b max ps 6.5 fm and E = 350 AMeV are considered. Fig. 4 illustrates this 
effect. The nucleons show positive p x - values for small rapidities on the average in 
the forward hemisphere (y cm > 0) whereas negative p x -values are observed for higher 
rapidities. This effect is seen for the hard equation of state without momentum de- 
pendent interactions, it is very sensitive to the incident energy, the impact parameter, 
and most importantly, to the addition of momentum dependent interactions (mdi). 
The signs of the average p^-values become positive for all positive rapidities if the im- 
pact parameter is reduced to b = 0.7b max . The same happens if momentum dependent 
interactions (which give additional repulsion) are introduced. The following scenario 
might explain the two components: Nucleons which have experienced higher densi- 
ties, e.g. p max > 1.3po are preferentially visible at small rapidities. This compressed, 
stopped matter shows positive flow. The spectator matter, which has experienced less 
compression, shows negative flow. The separation of the two components is clearly 
visible when applying a cut on the maximum density for slightly different system 
parameters (E = 330 AMeV and b = 0.75b max ). In addition the components can 
be separated with respect to their type in a simple configuration space coalescence 
model PH| . Protons yield the major part of the component at midrapidity whereas 
heavier fragments rule the outer component. The time-evolution of the collision can 
be imagined as if the spectators were sucked to the participant zone. 
A two-component flow is observed for Ca+Ca at 170 AMeV e.g. at b = 0.8b max with 
momentum dependent interactions, too. The sign of the components at midrapidity 
and larger rapidities is just opposite to those observed without mdi. Coalescence 
considerations indicate in turn that the components around midrapidity stem from 
heavier fragments, while free nucleons contribute mainly at y/y pro j. > ±1- 
This double sign change is highly sensitive to the momentum dependent interactions 
and should therefore be experimentally scrutinized. 
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Let us now turn to another point: A smaller balance energy E^i is expected for the 
heavy system Au+Au (A = 394) than for Ca+Ca due to the cited A~^-law. Experi- 



mentally so far only an upper bound for the balance energy of E^i < 60 AMeV fl2 
has been found. Therefore the existence of negative flow is an open question due to 
the strong Coulomb repulsion. We show that this is due to an ill-defined frame of 
reference. The flow is in fact balanced at E^i = 55 ± 5 AMeV and E^i = 65 ± 5 
AMeV for the impact parameters b = 0.256 max ~ 3.3 fm and b = 0.5b max w 6.5 fm, 
respectively and for a hard equation of state without momentum dependent inter- 
actions. These values are obtained if the initial pre-contact rotation of the system 
due to Rutherford-trajectories is subtracted. In this system the sign-reversal for the 
reduced flow is clearly visible. Fig. 5 shows the respective calculation for Au+Au at 
50 AMeV and b = 0.5b max ~ 6.5 fm. In the rotated system the flow is obviously 
negative whereas a flat distribution is obtained in the laboratory frame. In the lab- 
oratory frame negative flow does not appear for any impact parameter, even not for 
low energies. 

Let us now turn to the squeeze-out which is an established effect [|57|-|3"9"||. Excited 
participant matter is pushed out perpendicular to the reaction plane. At energies 
dicussed in this paper this behaviour might be different. In fig. 6 these azimuthal 
angular-distributions are plotted for the system Ca+Ca (hard eos+mdi) at their 
respective balance energies with different impact parameters. The considered ra- 
pidity is —0.15 < y/y p < 0.15 according to recent experiments for the heavier sys- 
tem Zn+Ni [jiOfl . The full lines are the result of fits by the Legendre-expansion: 
dN/dcp = a (l + aicos(0) + a 2 cos(20)). The value of a 2 gives a measure of the 
anisotropy of this collective motion. Negative values of a 2 show prefered emission 
perpendicular to the reaction plane whereas positive values describe an enhancement 
in the reaction-plane. Fig.6 shows that for Ca+Ca the in-plane emission is prefered 
for larger impact parameters, and a slight out-of-plane enhancement is observed for 
rather central collisions at the balance energies and at midrapidity. The transition 
energy where the anisotropy parameter a 2 becomes zero, corresponding to an az- 
imuthally symmetrical distribution, was measured for Zn+Ni . It was found that 
this transition energy is smaller than the corresponding balance energy. Our calcu- 
lations for the lighter system Ca+Ca show the transition energies to be larger than 
the balance energy for larger impact parameter (b > 0Ab max ), but smaller for more 
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central collisions. This was already indicated by measurements for Ar+V fH ]. Mea- 
surements indicate that the in-plane enhancement increases with impact parameter 



42fl . This can be seen in fig. 7 for Ca+Ca at 80 AMeV and various impact parame- 



ters. Light fragments show a slightly more pronounced in-plane to out-of-plane ratio 
than single nucleons if clustering is taken into account. Consequently, it must be 
pointed out that even at the in-plane balance energy collective flow characteristics 
are clearly visible in the azimuthal distributions. 



IV. CONCLUSIONS 



We have investigated the disappearance of the in-plane flow for Ca+Ca and 
Au+Au. 

•A strong impact parameter dependence of the in-plane balance energy is ob- 
served. The balance energy clearly increases with impact parameter. This cannot be 
neglected while pinning down basic properties of excited nuclear matter. 
•The balance energy is smaller with momentum dependent interactions than without 
for large impact parameters. The difference might be a tool to get information about 
the proper parametrization of the momentum dependent interactions. 
•Negative flow angles will not be visible in the laboratory frame for the heavy Au+Au 
system due to the long range Coulomb forces, although the in-plane flow disappears. 
Negative flow and the respective balance energies are visible in the frame where the 
pre-contact rotation due to the initial Rutherford-trajectories is subtracted. How- 
ever, a maximum mass must exist where negative flow can still be observed in the 
laboratory frame. 

•A new two-component flow was shown for large impact parameters. One com- 
ponent stems from participant particles at rapidities around y cm whereas the other 
component results from cold spectator matter. They show opposite sign in the p x {y)- 
distribution. The existence of two distinctly different flow-components depends on 
the inclusion of momentum dependent interactions. This is of great importance for 
the proper determination of the parametrization of the momentum dependent inter- 
actions or other basic properties such as the in-medium NN-cross section. 
•Finally, azimuthal distributions demonstrate the existence of flow, even at the bal- 
ance energy. For the system Ca+Ca the energy of the change from an preferentially 
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in-plane to out-of-plane emission is smaller for central collisions and larger for in- 
creasing impact parameters than the balance energy. This energy of an azimuthally 
symmetrical distribution can provide valuable information complementary to the 
in-plane balance energy. The in-plane to out-of-plane ratio increases with impact 
parameter. 

The search for tools to describe excited nuclear matter in nucleus-nucleus collisions 
and the search for signals to determine unambigiously the basic physical attributes 
is going on. 
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FIG. 1. Transverse momentum projected onto the reaction plane function of 

the normalized rapidity. This p x (?//yp)-distribution of protons for the system Ca+Ca is 
plotted at the two incident energies, 80 AMeV and 130 AMeV. The impact parameter is 
half the maximum impact parameter b = 0.5b ma x- For each curve thousand events were 
calculated with a hard equation of state without momentum dependent interactions. The 
lines are plotted to guide the eye. 
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FIG. 2. Reduced flow values as a function of incident energy and impact parameter 
for Ca+Ca. The impact parameters are 0.25, 0.4, 0.5, 0.6 times the maximum impact pa- 
rameter. Each point is a result of thousand events with a hard equation of state without 
momentum dependence. The straight lines are the results of linear fits. 
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FIG. 3. The in-plane balance energy Eb a i as a function of impact parameter b for 
the system Ca+Ca. The circels and squares are the calculated values without and with 
momentum dependent interactions, respectively. The curves are plotted to guide the eye. 




FIG. 4. p x (y/y p )-distribution of protons for the semi-peripheral (b = 0.85b max ) colli- 
sion of Ca+Ca at 350 AMeV incident energy. This two-component flow is received by a 
calculation of 10000 events with a hard equation of state without momentum dependent 
interactions. 
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FIG. 5. pa;(y/yp)-distribution of protons for the system Au+Au at 50 AMeV. The 
impact parameter is 0.5b max . The squares and circles correspond to calculations with and 
without an initialization on Rutherford trajectories. 500 events were calculated for each 
curve with a hard equation of state without momentum dependence. 
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FIG. 6. Azimuthal distributions with respect to the reaction plane for Ca+Ca. 
The incident energies and impact parameters correspond to the determined in-plane 
balance energies Eb a i(b) with momentum dependent interactions. The rapidity- 
range is restricted to —0.15 < y/y p < 0.15. The curves are fits according to 
dN/d(f) = ao(l + aicos((p) + a,2Cos(2(fi)). 
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FIG. 7. Azimuthal distributions with respect to the reaction plane for Ca+Ca at 80 
AMeV and for three different impact parameters b = 0.25,0.4, and 0.5b max . 



